Extreme tunability of interactions in a 7 Li Bose-Einstein condensate 
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We use a Feshbach resonance to tune the scattering length a of a Bose-Einstein condensate 
of 7 Li in the \F = l,mj? = 1) state. Using the spatial extent of the trapped condensate we 
extract a over a range spanning 7 decades from small attractive interactions to extremely strong 
repulsive interactions. The shallow zero-crossing in the wing of the Feshbach resonance enables the 
determination of a as small as 0.01 Bohr radii. Evidence of the weak anisotropic magnetic dipole 
interaction is obtained by comparison with different trap geometries for small a. 
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The ability to control the parameters of ultracold 
atomic gases and to impose external potentials upon 
them provides unique opportunities to create model sys- 
tems for exploring complex phenomena in condensed 
matter and nuclear physics. Control of atomic interac- 
tions using Feshbach resonances has proven to be partic- 
ularly productive in studies involving Bose-Einstein con- 
densates (BECs) or paired Fermi gases While strong 
interactions are usually the focus of these studies, inter- 
esting phenomena also occur in the weakly-interacting 
regime. An example of such a phenomenon is Anderson 
localization in disordered media [2|, which was recently 
observed in weakly repulsive BECs 0, 0| ■ Another ex- 
ample is the formation of bright solitons in BECs with 
weakly attractive interactions, which have been created 
in condensates of 7 Li 0, [f| and 85 Rb Q. Atom inter- 
ferometers may also benefit by the increased coherence 
times afforded by weakly interacting gases IRQ, or even 
by a dispersionless atomic soliton laser 0, [lfj ] . 

Several atomic species exhibit Feshbach resonances 
where the s-wave scattering length a changes sign at a 
certain field in the wings of the resonance. These zero- 
crossings are useful in applications requiring weak inter- 
actions. In addition to 7 Li [1,0], such zero-crossings have 



been studied in 85 Rb [lj, 52 Cr [I2j, 39 K [Pj, and 133 Cs 
Q, In this Letter, we report the measurement of a for 
7 Li in the \F = 1, mj? = 1) state for fields near the Fes- 
hbach resonance at 737 G 0, @, [3, 15 1. By measuring 
the in situ size of the confined condensate, a is mea- 
sured over a range of 7 decades. We find that the slope 
of the zero-crossing is only ~0.1ao/G, where ao is the 
Bohr radius. This is the shallowest known zero-crossing, 
requiring only modest field stability to achieve an essen- 
tially non-interacting gas. We explore the effects of the 
magnetic dipole interaction (MDI) in this regime. Unlike 
Cr, which has a large magnetic moment of 6 resulting 



in a relatively large MDI 0, [TJ [3 , the MDI in alkali 
atoms is weak due to their small magnetic moments of 
~1 /ljb- Nonetheless, the MDI has been recently detected 
in alkali atoms IjJ 2(| ■ We explore the role of the MDI 
by modifying the confining geometry of the BEC. 

Our experimental apparatus for generating a BEC of 
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7 Li has been described previously [6|, l2l| . Atoms in the 
\F = 1 , tti p = 1) state are confined in an optical trap 
formed from a single focused laser beam with wavelength 
of 1.03 ^m. A bias magnetic field, directed along the trap 
axis, is used to tune a via the Feshbach resonance. We 
create condensates at a field where a is large to facili- 
tate rapid rethermalization of the atoms during evapo- 
ration from the optical dipole trap. After a condensate 
is formed we slowly (~4s) ramp the field to the desired 
value and determine the scattering length, as described 
below. There is no discernable thermal part to the den- 
sity distributions and we estimate that T/T c < 0.5, where 
T c is the condensation temperature. The final trap- 
ping potential is a combination of the optical field and 
a residual axial magnetic curvature from the bias field. 
The trap is cylindrically symmetric with measured ra- 
dial and axial trapping frequencies of uj r /2ir = 193 Hz 
and w 2 /27r = 3 Hz, respectively. 

We use in situ polarization phase-contrast imaging 22 1 
to acquire the column density distribution of the conden- 
sate at the desired magnetic field. When the s-wave in- 
teractions are large and repulsive they inflate the size of 
the condensate well above the harmonic oscillator size. 
As the interactions decrease the size of the condensate 
becomes smaller, approaching the harmonic oscillator 
ground state near zero interactions. Figure [T] shows rep- 
resentative images of condensates with various repulsive 
or attractive interaction strengths. Solitons form when 
a < 0, either a single one for a slow magnetic field ramp 
or multiple solitons for ramps fast compared to the axial 
trap period. 

We integrate the image of the condensate in the re- 
maining radial dimension to produce an axial density 
profile. The 1/e radius of this profile is used as a measure 
of the condensate size, as shown in Fig. [2] for a range of 
magnetic field values. In the Thomas-Fermi regime, the 
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FIG. 1: (color online) Representative in situ polarization 
phase-contrast images of condensates with various interaction 
strengths, (a) B = 719.1 G, a = 396a , N = 1.7 X 10 5 ; (b) 
B = 597.4 G, a = 8a , N = 2.9 x 10 5 ; (c) B = 544.7 G, a = 
0.1 oo, N = 2.0xl0 5 ; (d) B = 542.4 G, a = -0.1 oo, N = 1.2x 
10 5 ; (e) same as (d) but with a faster field ramp from 710 G 
to 542.4 G, resulting in multiple solitons with N « 10 4 per 
soliton. The probe laser detuning from resonance is adjusted 
to keep a nearly constant signal level, and varies between 
20 7 for large a to 150 7 for small a, where 7/2-7T » 5.9 MHz 
is the excited state linewidth. The color map is adjusted to 
maximize contrast for each image. 

axial size of the condensate is dependent on the product 
of a and the number of atoms in the condensate N. The 
average number per condensate is Nq — 3 x 10 5 atoms, 
with a shot-to-shot variation of 20%. The inset of Fig. fj] 
shows the axial size scaled by (N/Nq) 1 ^ 5 to account for 
these fluctuations. Several condensates are found to have 
axial sizes smaller than the axial harmonic oscillator size 
due to net attractive interactions, as discussed below. 

To determine a for each image requires a mapping from 
the measured axial size and N to a. We model the sys- 
tem using the three-dimensional (3D) Gross-Pitaevskii 
equation 

ft 2 ^9 , t, , 4:TtH 2 a . , .0 , 

fj,ip= - — V 2 ip + Vijj + \ip\ 2 ip 

2m m 

where we account for the MDI in addition to the s-wave 
contact interaction and the trapping potential. At 540 G, 
fx m w 0.94 fi B for 7 Li in the |1, 1) state. Mapping is ac- 
complished by performing a variational calculation using 
a 3D cylindrically symmetric Gaussian wavefunction as 
a trial solution to Eq. |T]). Minimizing the correspond- 
ing energy functional results in equations for the radial 
and axial sizes of the condensate [23j , which are solved to 
give the desired mapping function. Figure [3] shows this 
mapping with and without inclusion of the MDI, as well 
as the corresponding mapping from the Thomas-Fermi 
approximation. In our geometry, the magnetic moments 
are aligned with the long axis of the trap. This causes the 
MDI to be effectively attractive, making the condensate 
smaller axially for a given value of a. We have verified the 
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FIG. 2: (color online) Axial size of the condensate as a func- 
tion of magnetic field. The axial size is defined as the 1/e 
radius of the axial density profile and is scaled by the axial 
harmonic oscillator size l z — y/fi/mtxiz « 22 /im. The resolu- 
tion of the optical imaging system is ~3.3 /im (dotted line). 
The dashed line is the size of the condensate (I4 « 0.62 l z ) 
found by solving Eq. ([l]) with a = 0. The zero-crossing in a 
occurs when the size of the condensate equals Id (arrow). Ne- 
glecting dipolar effects results in a zero-crossing about 0.5 G 
higher, where the axial size equals Z z . The inset shows the 
axial size corrected for number variation as described in the 
text. Individual data points and error bars are the average 
and standard error of approximately 10 shots taken at each 
field. Systematic uncertainty in the axial size is ~3% from 
uncertainty in temperature and the uncertainty in imaging 
magnification. The systematic uncertainty in magnetic field 
due to calibration (via radio frequency transitions from the 
|2, 2) to the |1, 1) state) is ~0.1 G. We have binned the data 
into intervals of this size. 



accuracy of the variational calculation by exact numerical 
solution of Eq. ([T]) for various values of a and find good 
agreement between the two methods. Since the varia- 
tional calculation is much faster computationally, we use 
this method to analyze the data. 

Figure [4] shows the axial size data of Fig. [2] mapped 
onto a. The general shape follows that of a typical 
Feshbach resonance with a = obg[1 + A/(£? — Boo)], 
where a BG = -24.5^ a , A = 192.3(3) G, and = 
736.8(2) G. The uncertainties in these derived values 
are a result of the systematic uncertainty in field cal- 
ibration of 0.1 G and a systematic uncertainty in a of 
^20%, primarily due to uncertainty in measuring the 
axial size and determination of u> z . A linear fit to the 
data for B < 550 G gives a slope of 0.08(1) ao/G and 
a zero-crossing at Bq — 543.6(1) G [24|. The smallest 
mean positive scattering length of a collection of shots 
was 0.01(2) a at 543.6(1) G with -3 x 10 5 atoms. Under 
these conditions the peak density is 3 x 10 14 cm -3 and the 
corresponding condensate healing length is comparable 
to the length of the condensate itself. Although Eq. ((T|) 
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FIG. 3: (color online) Mapping functions of axial size to 
a using a Gaussian trial wavefunction in a variational solu- 
tion to Eq. {l}, including (solid) and neglecting (dashed) the 
MDI; also shown is the Thomas- Fermi approximation (dot- 
ted). These mappings were computed for N = 3 x 10 5 , 
uj r /2TT = 193 Hz and u> z /2n — 3 Hz. In practice, we compute 
the mapping individually for each imaged condensate to ac- 
count for variations in N and a field dependent variation in ui z 
of ~5% over the relevant magnetic field range. The Gaussian 
solution neglecting the MDI asymptotically approaches l z at 
zero interactions, while their inclusion causes the solution to 
asymptotically approach a value smaller than l z . 



assumes the mean field approximation, beyond mean field 
corrections are expected to be important when na 3 > 1 
H, [H, [27|, [H, Est The leading order correction to the 
interaction term in Eq. ([TJ) , the Lee-Huang- Yang param- 



eter, is a = 32/(3 v / 7r)" i 



^> 1 for the most strongly 



interacting condensates observed. We have accounted 
for this correction in extracting a for data where a < 1. 
For the four data points with a > 1, this low-density 
expansion is not valid. We are unaware of a theoretical 
treatment that addresses the density distribution in the 
strongly interacting regime. While we extract a value a 
a for these four data points by fitting to a Thomas-Fermi 
profile ignoring beyond mean-field effects, and plot them 
in Fig. 4, we exclude them in the Feshbach resonance 
fit. Using this method the largest mean positive scatter- 
ing length was -2 x 10 5 a at 736.9(1) G with -2 x 10 4 
atoms, which has a peak density n « 5 x 10 10 cm~ 3 . 
The comparatively smaller number of atoms close to res- 
onance is likely due to large inelastic collisional losses in 
this regime [30|. 

Figure 2] also shows a comparison between a coupled- 
channels calculation and the experimentally extracted 
values of a. The coupled-channels calculation requires 
the ground-state singlet and triplet potentials of 7 Li2 
as input, as described previously [3lL 132^. We have up- 
dated the potentials to be consistent with the new mea- 
surements of and Bq reported here, as well as the 
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FIG. 4: (color online) Axial size data of Fig. [2] mapped onto 
a. Results of a coupled-channels calculation are shown by 
the solid line. The Feshbach resonance fit is indicated by the 
dashed line. The inset shows the extracted values of a near the 
zero-crossing. The mean and standard error of approximately 
10 shots taken at each field is shown. In addition, we estimate 
a systematic uncertainty of ~20% in a. 



previously measured bin ding energy of the least bound 



triplet vibrational level [31|, |33j. The updates involve 
adjusting the singlet and triplet dissociation energies 
to D e (X 1 S+) = 8516.68(10) cm' 1 and L> e (a 3 £+) = 
333.714(40) cm" 1 , where the stated uncertainties account 
for uncertainties in the remaining portions of the model 
potentials. These values are consistent with previous de- 
terminations 32j, [3J, [35|] . The agreement between the 
calculated and measured values of a, while not perfect 
over the entire range of fields, is reasonably accurate over 
a range spanning 7 decades. 

The effects of the MDI are strongly dependent on ge- 
ometry. To better distinguish their role, we increased 
the axial trapping frequency from 3 Hz to 16 Hz by ap- 
plying magnetic curvature. Figure [5] compares the ex- 
tracted values of a for both trap geometries when the 
MDI is included or neglected in the mapping function. 
As expected, neglecting the MDI in the analysis system- 
atically lowers the extracted values of a. This effect is 
most noticeable near the zero-crossing where a system- 
atic geometry-dependent discrepancy appears in the de- 
rived values of a. Inclusion of the MDI in the analysis 
produces a consistent value of a for a given magnetic 
field regardless of the trapping potential. The data show 
that the magnetic dipole interaction, although weak, is 
discernible in 7 Li despite having a magnetic moment of 
only ~l/x B . 

We have mapped the Feshbach resonance from the 
regime of small attractive interactions far from the res- 
onance to extremely strong repulsive interactions very 
close to resonance. The zero-crossing and resonance posi- 
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FIG. 5: (color online) Extracted values of a near the zero- 
crossing for trapping potentials with u} z /2n = 3 Hz (filled 
squares) or u) z /2ir = 16 Hz (unfilled squares), when (a) ne- 
glecting or (b) including the MDI in the mapping function. 
The MDI has a negligible effect on the extracted values of 
a for the 16 Hz trap, but neglecting the MDI in analysis of 
the 3 Hz trap systematically lowers the mapped values of a, 
especially for a < 0.15 ao. 



tions have been precisely located, enabling experimental 
access to a broad range of accurately known interactions. 
Of particular interest will be explorations of atom and 
soliton transport through a disordered potential in the 
weakly interacting regime. 
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